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We study the glued-trees problem of Childs, et. al. [T] in the adiabatic model of quantum com- 
puting and provide an annealing schedule to solve an oracular problem exponentially faster than 
classically possible. The Hamiltonians involved in the quantum annealing do not suffer from the so- 
called sign problem. Unlike the typical scenario, our schedule is efficient even though the minimum 
energy gap of the Hamiltonians is exponentially small in the problem size. We discuss generaliza- 
tions based on initial-state randomization to avoid some slowdowns in adiabatic quantum computing 
due to small gaps. 
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Quantum annealing is a powerful heuristic to solve 
problems in optimization [21 [3] . In quantum computing, 
the method consists of preparing a low-energy or ground 
state of a quantum system such that, after a sim- 
ple measurement, the optimal solution is obtained with 
large probability. is prepared by following a particu- 
lar annealing schedule, with a parametrized Hamiltonian 
path subject to initial and final conditions. A ground 
state of the initial Hamiltonian is then transformed to \ip) 
by varying the parameter adiabatically. In contrast to 
more general quantum adiabatic state transformations, 
the Hamiltonians along the path in quantum annealing 
are termed stoquastic and do not suffer from the so-called 
numerical sign problem [2J: for a specified basis, the off- 
diagonal Hamiltonian-matrix entries are nonpositive [5]. 
This property is useful for classical simulations [3]. 

A sufficient condition for convergence of the quantum 
method is given by the quantum adiabatic approxima- 
tion. It asserts that, if the rate of change of the Hamil- 
tonian scales with the energy gap A between their two 
lowest-energy states, can be prepared with controlled 
accuracy [5J [7] . Such an approximation may also be nec- 
essary [5J. However, it could result in undesired over- 
heads if A is small but transitions between the lowest- 
energy states are forbidden due to selection rules, or if 
transitions between lowest-energy states can be exploited 
to prepare The latter case corresponds to the anneal- 
ing schedule in this Letter. It turns out that the relevant 
energy gap for the adiabatic approximation in these cases 
is not A and can be much bigger. 

Because of the properties of the Hamiltonians, the an- 
nealing can also be simulated using probabilistic classical 
methods such as quantum Monte-Carlo (QMC) [9]. The 
goal in QMC is to sample according to the distribution 
of the ground state, i.e. with probabilities coming from 
amplitudes squared. While we lack of necessary condi- 
tions that guarantee convergence, the power of QMC is 
widely recognized [3J [21 [TO] - In fact, if the Hamiltonians 
satisfy an additional frustration-free property, efficient 
QMC simulations for quantum annealing exist [111 1 1 2] . 
This places a doubt on whether a quantum-computer 
simulation of general quantum annealing processes can 



ever be done using substantially less resources than QMC 
or any other classical simulation. 

Towards answering this question, we provide an oracu- 
lar problem and give a quantum- annealing schedule that, 
on a quantum computer, prepares a quantum state \tp) 
encoding the solution. The time required to prepare \tp) 
is polynomial in the problem size, herein poly(n). The 
oracular problem was first introduced in Ref. [1] in the 
context of quantum walks, where it was also shown that 
no classical method can give the solution using poly(n) 
number of oracle calls. Our result thus places limits on 
the power of classical methods that simulate quantum 
annealing, even when the sign-problem is not present. 

We remark that the general question of existence of 
efficient classical simulations when A is l/poly(n) is not 
answered in this Letter. The annealing schedule we pro- 
vide is not intended to follow the ground state in the 
path; transitions to the closest (first-excited) eigenstate 
are allowed. Nevertheless, the system (almost) remains 
in the subspace spanned by these two states at all times. 
There are regions in the path where A oc exp(— n). We 
induce transitions in that subspace by choosing an an- 
nealing rate that is much larger than A, i.e. at l/poly(n) 
rates. Contrary to the typical case, such transitions are 
useful here. They guarantee that I?/ 1 ) is prepared after 
the annealing due to a symmetry argument: The same 
type of transition that transforms the ground to the first- 
excited state, later transforms the first-excited state back 
to the final ground state 

In more detail, we consider the oracular problem from 
Ref. [T] that is defined as follows. We are given an oracle 
that consists of the adjacency matrix A of two binary 
trees that are randomly glued (by a random cycle) as in 
Fig[l] Specifically, there are N € 0(2") vertices named 
with randomly chosen 2n-bit strings. The oracle outputs 
the names of the adjacent vertices on any given input ver- 
tex name. There are two special vertices, ENTRANCE 
and EXIT - the roots of the binary trees. They can be 
identified because they are the only vertices of degree two 
in the graph. The problem is: Given an oracle A for the 
graph and the name of the ENTRANCE, find the name 
of the EXIT. As mentioned, no classical algorithm can 
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FIG. 1: Two binary trees of depth n = 4 glued randomly. The 
number of vertices is TV = 2 n+2 — 2. Each vertex is labeled 
with a randomly chosen 2n-bit string, j is the column number. 



output the name of the EXIT, with high probability, us- 
ing less than a subexponential (in n) number of oracles. 
A quantum walk-based algorithm can solve this prob- 
lem efficiently with bounded probability [TJ. However, 
a quantum annealing method to solve this problem effi- 
ciently using simple (stoquastic) Hamiltonians remained 
unknown. In this Letter, we present a new quantum an- 
nealing approach for this problem that provides the so- 
lution with arbitrarily high probability, c.f. Refs. [71 113j. 

We assume a Hamiltonian version of the oracle so that 
evolutions under A can be implemented. We also allow 
for evolutions under H and Hi, these being the diago- 
nal Hamiltonians that distinguish the ENTRANCE and 
EXIT, respectively. Such evolutions can be realized ef- 
ficiently [21 [T5], i.e. using G(T) oracles for evolution 
time T > 0. We let a(V) G {0, l} 2 ™ be the name of 
vertex V. Then, H \a') = -^(entrance), a' W) and 
Hi \a!) = — 5 a (ExiT),a' l a ')) so that their ground states 
encode a(ENTRANCE) and a(EXIT), respectively. The 
Hamiltonian path for the annealing will consist of a spe- 
cific interpolation involving Hq, A, and Hi. 

As in Ref . pQ , we find it useful to define the (orthonor- 
mal) states 



| col,-) = 



E 



*(*)> 



(i) 



^ i£jth column 



These are uniform-superposition states over all states la- 
beled by the names of vertices at the j— th column. Nj — 
23 for < j < n and Nj = 2 2n+1 ^ for n+1 <j< 2n+l; 
see Fig. IT) In particular, |col ) = | a (ENTRANCE)) and 
|col2„+i; = | a (EXIT)). We observe that the subspace 
spanned by {|colj)}o<j<2n+i is invariant under the action 
of A, H , and Hi. In the basis determined by Eqs. ()TJ), 
A has non-zero matrix elements in its first off-diagonals 
only. For simplicity, we redefine A <— \[2A so that the 
matrix elements arc 



(colj| A|col i+ i) = 



V2 j — n 
1 otherwise. 



(2) 



Also, 

(coljl H Q \co\j) = -S jfi , (3) 

(colj| Hi \colj) = Sj^n+l ■ 

We choose the Hamiltonian path 

H(s) = (l-s)aH -s(l-s)A + saHi (4) 

that interpolates between Ho and Hi for < s < 1. 
The parameter a is independent of n and satisfies < 
a < 1/2. We will show that using the Hamiltonian path 
of Eq. Q, and annealing at a rate s(t) cx l/poly(n), 
the resulting evolution transforms | a (ENTRANCE)) to 
a state that has arbitrarily high overlap with | a (EXIT)). 
Spectral properties — To prove the efficiency of the quan- 
tum method we utilize the spectral properties of H(s); 
particularly relevant arc the spectral gaps. The following 
analysis is valid if we restrict to the invariant subspace 
spanned by Eqs. (JlJ). Figure [2] shows the three lowest 
eigenvalues of H(s), obtained numerically, in this sub- 
space. This suggests a particular eigenvalue behavior. 
We can analytically study the Hamiltonians by propos- 
ing the ansatz \<p) = Y]. jj |colj), with 



ae 



ipj 



+ be~ lpl , < j < n , 



(5) 



ce 



ip(2n+l-j) + de -ip&n+l-j) , n + 1 < j < 2n + 1 , 



and p € C. The eigenvalue condition H\<j>) — \ \4>) and 
(cj}\ <fi) = 1 allow us to find expressions for a, b, c, d, and 
A. In particular, the eigenvalues are A = — 2s(l — s) cosp. 
We provide a more detailed analysis of the spectrum in 
the supplemental online material (Supp. Mat.) and 
present only the relevant results here. Because of the 
so (1 — s) symmetry, it suffices to analyze the parame- 
ter region < s < 1/2. 

In the following, x « e y if \x — y\ < e and e G 0(2~™/ 2 ). 
When n — y oo, the two lowest eigenvalues cross, or be- 
come equal, at s — s x = aj\j2. Different eigenvalue 
behavior is obtained at both sides of s x . For n < oo 
and < s < s x , the spectral gap between the two lowest 
eigenvalues is 



A 10 (s) = Ai(s) - A (s) 
3s 



-(1 



(G) 



The eigenvalue crossing is avoided for n < oo but the 
spectral gap Aio(s) is exponentially small in n near s x - 
Also, Aio(s) > c/n 3, for < s < Si = s x — 5, with 
S G fi(l/n 3 ). Because the ordering of the two lowest 
eigenvalues swaps for s > s x , Aio(s) > c/n 3 in the region 
s 2 = s x + 6 < s < 1/2. 

The second excited state has an eigenvalue that corre- 
sponds to p m v n/(n + 1), with v G 0(l/(n + l) 2 ). The 
spectral gap with the first excited state for < s < a is 

A 21 (s) = A 2 (s)-A 1 (s) 

s(l - s) (cos (_J T ) - A j . (7) 
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FIG. 2: The three lowest eigenvalues of H(s) in the subspace 
spanned by the states of Eqs. (JTJ, for a = l/y/8 and n = 10. 
Aj'fc — ~ is the gap between Xj and Afc, the eigenvalues 
corresponding to the j-th and fc-th excited states, respectively. 
We divide the evolution in five stages according to si, S2, S3, 
and S4 (see text), with si < s x = a/y2 = 0.25 < S2 and 
S3 = l — S2<1 — s x < 84 = 1 — si- Inside [81,82] and [S3, 84], 
the gap Aio(s) becomes exponentially small in n. Elsewhere, 
A10 is only polynomially small in n. Brown arrows depict level 
transitions for an annealing rate in which s(t) oc l/poly(n). 
Other gap scalings are also shown. 

In particular, A 21 (s) G 0(1) in the region si < s < a. 

For a < s < 1/2, the second eigenvalue also corre- 
sponds to p w/(n + 1). To bound the gap with the 
third eigenvalue, a more detailed analysis that approxi- 
mates p at order l/(n + l) 2 is carried in the Supp. Mat.. 
It results in A 2 i(s) > c'/n 3 , for some c' > (see Fig. [2]). 
Annealing schedules — We use the following adiabatic 
approximation from Refs. [6]. Let the initial state be 
an eigenstate of H(s ) and A(s) the spectral gap to 
the nearest (non-degenerate) eigenstate in some region 
so < 9 < Sf. Then, an annealing rate of s(t) oc eA 2 (s) or 
smaller suffices to prepare the eigenstate of H(sf) at er- 
ror amplitude smaller than y/e(sf — sq); i.e. the overlap 
between the evolved state and the eigenstate is at least 
a/1 — e(sf — sq). (Better error scaling is possible [71 fT3"].) 

To prove that s(t) oc e/n e suffices to transform 
ja(ENTRANCE)) to |o(EXIT)) with large overlap, we 
split the evolution according to [0,1] = Ui=i^i' with 
Vi = [0, sx), V 2 = [si,s 2 ), V 3 = [82,83), V 4 = [s 3 ,s 4 ), 
and V5 — [s 4 , 1]. The values of were determined pre- 
viously; see Fig. [2} We write \4>o(s)) and \(f>i(s)) for the 
ground and first excited states, respectively. Note that 
|0 O (O)) = |a(ENTRANCE)) and |0 O (1)) = |o(EXIT)). 
Then, due to the gap bounds and the adiabatic approxi- 
mation, the following transformations occur: 

l<M0)) \M*i)) > (8) 

—} x indicates that the transformation occurred at error 



amplitude of order x. 

Because A 2 i(s) G f2(l) for s G V%, transformations 
between the ground or first-excited state and the sec- 
ond excited state occur with amplitude smaller than y/e. 
Thus, all relevant transitions in V2 occur in the manifold 
spanned by {\4>o(s)) , \4>i(s))}. For our annealing rate, 
the following transformations occur with large amplitude 
(see below): 

\M«i)) \<PM) , (9) 

I0l(s 3 )) -> \<t>o( s *)) ■ 

To obtain the approximation errors for Eqs. @, we 
introduce the state \u) that is a uniform superposition 
over all vertex names: 

2n+l rjT~ 

M = 7m E l«W)=E V^ |co1 ^ • (10) 

iGgraph j=0 

Here, Nj = "2? for < j < n and Nj = 2 2 "+ 1 "-' for 
n + 1 < j < 2n + 1. Interestingly, \u) is almost an 
eigenstate for all s: H(s) \u) « e / 2 — (s(l — s)3/v / 2) \u) 
and e G C(2-™/ 2 ) (see Supp. Mat.). We define f(t) = 
I (u| t7(i) \u) | 2 , where L7(i) is the evolution operator and 
/(0) = 1. Schrodinger's equation yields 

f(t) = -i (u\ H(s(t))U(t) \u) (u\ U\t) \u) + c.c. 

« e 0. (11) 

If T G 0(n 3 ) is the evolution time to change s from si 
to s 2 with our annealing schedule (|s 2 — si| € fi(l/n 3 )), 
we have /(T) sa e / 1 for e' G 0(en 3 ). In addition (Supp. 
Mat.), 

Ku|fc(ai))|««' 1, |H^ (s 2 ))|«e'l, (12) 
resulting in 

I (^i(si)| Z7(T) |^o(s2)) | 2 «5e' 1 ■ (13) 

The transformation |</>i(si)) — > y-gjy | <^o (• s 2 )) then oc- 
curs. Moreover, because U(t) is unitary, we also have 
\<Po( s i)) ~ ► 1^1(52))) and from symmetry arguments, 
I^i(s3)) — ^y5^_ |0o(-S4)) - These level transitions are 
shown in Fig. [2j Together with the transformations in 
Eqs. ([8]), they prove the success of our quantum anneal- 
ing method. Because e' <C 1 for large n, the overall 
amplitude error is dominated by that of the adiabatic 
approximation. This is of order y/e, with e arbitrary. 

In Fig. [3] we show the overlaps of U(t) |0q(O)) with 
the ground and first excited states as a function of the 
evolution time, using our s(t), showing evidence for the 
transition among the two low energy levels. 
Initial state randomization: a generalization — A method 
that guarantees a successful final state preparation with 
probability 1/2 in this case is possible if we random- 
ize the initial state preparation: randomly prepare 
ja(ENTRANCE)) or \u) that is almost the first excited 
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FIG. 3: Overlap of the evolved state U(t) |a(ENTRANCE)) 
with the ground and first excited states as a function of time. 
In this case, n = 40, a = l/VE, and s(t) = t/10000. The 
simulation was performed in the column subspace. 



state. Both states can be prepared efficiently. If we 
choose an annealing rate s(t) oc e/n 6 that forbids transi- 
tions to the second excited level, we are guaranteed, from 
unitarity, that the final state is prepared with probability 
1/2. This method takes advantage of the fact that the 
gaps with the third eigenvalue are of order 1/n 3 in the 
region of interest. 

This randomized method can be carried to those cases 
in which the k lowest eigenvalues are at gaps l/poly(n) 
with the remaining ones. Then, any eigenstate in this 
manifold can be prepared with probability of order 1/k 
by randomizing the initial state preparation and choos- 
ing an adequate annealing schedule. If k is not too big, 
the method may be efficient regardless of whether the k 
lowest eigenvalues have small gaps or cross. It assumes 
the efficient preparation of any of the k initial states. 

General Hamiltonians do not satisfy the assumptions; 



however, some Hamiltonians of physical systems sat- 
isfy them. An example is the one-dimensional spin- 
1/2 Ising model in a constant transverse field: H(J) — 
JJ2i <T l <T z +1 + a xi with o~ l a being the Pauli operator on 
spin i. As the Ising coupling J is changed from to 
\J\ 1, the two lowest eigenvalues have a spectral 
gap that decreases exponentially with the system size 
n. However, the third eigenvalue is always at a distance 
l/poly(n) from the two lowest ones. This is a common 
property for systems that have conformal invariance. 

Conclusion — We provided an example of an oracular 
problem for which a quantum adiabatic evolution can 
find the solution efficiently whereas exponential time is 
required for any classical method. The Hamiltonians in 
the evolution do not suffer from the sign problem. Our 
result is a step towards proving the power of quantum 
methods for quantum annealing. We discussed the de- 
tails of why our algorithm works and how it can be gen- 
eralized, under some assumptions on the spectrum, by 
randomizing the initial state preparation. 
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I. SUPPLEMENTARY ONLINE MATERIAL 



We analyze the spectrum of 

H(s) = a(l - s)H - s(l - s)A + asHx 



(14) 



for < s < 1 and a = ~ 0.35 for simplicity (see 

below). The matrix elements of A are in Eq. (pi). We 
rescale the Hamiltonian and study H'(s) = H(sj/(s(l — 
s)) = a'Ho-A+p'Hi, with a' = a/s and /?' = a/(l-s). 
We let \4> p (s)) be an eigenstate and propose the ansatz 
\M S )) = E?"o 1 7? Icolj) given by 



7j = ae w 



7j 



be~ im , < j < n , (15) 
- de- lp3 ' , n + 1 < j < 2n + 1 , 



with f = In + 1 — j. The constants a, b, c, d are deter- 
mined below, p € C so that different solution behavior is 
obtained for real or imaginary p. 

From the eigenvalue equation H'(s) \4> p (s)) = 
A'(s) \4> p (s)), we obtain 

A' 

and the conditions 



-2cosp, (16) 



a (a + b) 
P'{c + d) 
ae *P(«+l) + he ~ip(n+l) 
ce v(n+l) + de -ip(n+l) 



ae~ lp + be ip , (17) 
ce- ip + de lp , 

V2(ce ipn +de- tpn ) , (18) 
V2(ae ipn + be~ lpn ) . 



Together with the normalization condition 1 = (</>|0), 
these determine the unknown parameters p,X,a,b,c,d. 
Simple inspection of Eqs. (17 1,(18) implies that p must 



obey the quantization condition 



f(p,n,a')f(p,n,(3 r )=2, 



with 



s'm((n + 2)p) — 9 sin((n + l)p) 
sin((n + l)p) — 9 sin(np) 



(19) 



(20) 



The symmetry a 1 f-> f3' is evident in Eq. ( 19 ) and thus 



it suffices to consider a',/3' for which < s < 1/2. The 
nature of the solutions depends on whether a > 1/2 or 
a < 1/2. Our choice of a = 1/V8 < 1/2 will produce 
an efficient annealing method. In the following, we use 
I to understand the behavior of 
ues and to estimate the scaling 



the condition of Eq. (19 



the three lowest eigenva 
of relevant energy gaps. 



1. Hyperbolic solutions 



We consider first the case where p = —iq is purely 
imaginary, so that ship — > sinhg and cosp coshg in 



G 



Eq. (pOb. We write f h (q,n,d) = f(-iq,n,6). The be- 



havior of the lowest eigenvalues for n < oo will be un- 
derstood by considering the large-n limit. In particular, 
lirrin^oo sinh((?i-|-l)(7o) — c/ sinh(ngo) = for <jo = ln(a'). 
Also, lim„_ j . oo // l (g ,w,/9') = a' if s < a < 1/2 (or 
a' > 1), which implies j3' < a'. Both limiting Eqs. im- 
ply that q determines an eigenvalue in the large-n limit. 
Another eigenvalue can be obtained by noticing that 



\u) is the uniform superposition state over all vertex 
names, i.e. 



2n+l 



\u) 



"iG graph 



3=0 



E i«w)-E\/^^i co1 .) > ( 27 ) 



lim fh{qi,n,a')fh(qi,n,P') 



(21) 



with Nj = 2 j for < j < n and Nj = 2 2 ™ +1 ^ for 
n+ 1 < j < 2n+ 1. Finite size corrections to \<fiq(s)) can 
be obtained from Eqs. (17) and (18). By simple matrix 



if qi = In \[2. The eigenvalues of the original Hamilto- 
nian of Eq. dT4l) are (s < 1/2) 



F(s) 
G(s) 



-2s(l — s) cosh qo 
-2s (1 — s) coshg! 



-{l-s){a + s 2 /a) 
-s(l - s)3/V2 , 



multiplication we obtain (a < 1/2) 

H(s) \u) « e/a G(s) \u) , (28) 
showing that \u) is almost an eigenstate for all s. Also, 



(«|fc(o))| = K«|0 ff (i))|« 2e i 



(29) 



respectively. Note that F(s) = G(s) if s = s x = a/\/2, 
where a level crossing for n — > oo occurs. 

The eigenvalues for n < oo are slightly different from 
those for n < oo. We are particularly interested in the 
energy gap around s x . For n — oo, a solution to Eq. ( 19 ) 
is obtained if we let q = In \[2 + e(s) and do a second- 
order Taylor series approximation; i.e. we approximate 



which follows from the spectrum of Hq and H\ , and the 
orthogonality condition. For other values of s, there is 
an energy gap A(s) between \<pq(s)) and any other eigen- 
state. Recall that the eigenvalue is ~ t /2 G(s). With no 
loss of generality \u) = o\<j)q(s)) + \/l — o 2 |_L), where 
(_L| <fiq(s)) — 0. We rescale the Hamiltonian by — G(s) 
so that H(s) \u) sa e / 2 and H(s) \4>q(s)) 



e/2 



2sinh(nq) « (V2) n {l+ne{s)_Hi 2 e 2 (s)/2-2 "). Plugging implies that yr^2 A ( s ) « £ or, equivalent^, 



the approximation in Eq. (|19j), we obtain two solutions 
with e(s) < e € C(2~ n / 2 ). This finite size correction 
also carries to the eigenvalues, proving an exponentially 
small gap of order e near s x . Away from s x , finite size 
corrections can be shown to be of smaller order, if s < a, 
by using a first-order Taylor series approximation. 

We write x ~ e y if | x — y\ < e. Then, for s < s x , the 
lowest eigenvalue for n < oo is 



0. This 



(30) 



Ao(s) w e/2 = -(1 - s){a + s 2 /a) 



(22) 



Similarly, the eigenvalue for the first excited state in that 
region is 



Ai(s) w e/2 = -s(l - s)3/V5 



(23) 



^. Gomometric solutions 

For a full characterization of the three lowest eigenval- 
ues, we investigate the second excited state for < s < 
1/2 as well as the first excited state for a < s < 1/2. 
To this end, we consider the eigenvectors of the gonio- 
metric type, i.e. those for which p S R, and consider 
the large-n limit first. Equation (20) is singular for 
P ^pr- Finite size corrections can be obtained by 



These two eigenvalues are swapped for s > s x (a level 
crossing in n — > oo). It implies that the lowest eigenvalue 
for 1/2 > s > s x is 



Ao(s) « e /2 G(s) = -s(l - 5)3/^2 . (24) 
The eigenvalue for the first excited state in a > s > s x 



allowing p(s) = + and using Eq. ( 20 ) to solve 

for x(s). By inspection, solutions for k = 1 correspond 
to the first and second excited eigenstates. We are ulti- 
mately interested in spectral gaps as they guarantee the 
convergence of the quantum annealing method. If A2(s) 
is the eigenvalue of the second excited state, it follows 
that, for < s < a, 



A 2 l(s) « e /2+i, s(l - s) 



Ai(s) w e/2 F(s) = -(1 - s)(a + s 2 /a) 



(25) 



-2 cos 



1 



3 

7^ 



The solution corresponding to q = \ny/2 + e(s) is valid 
for all s, i.e. it is the first excited state for s < s x and 



the ground state for s > s x . Using Eqs. (n7j) and (18), 



the corresponding eigenstate |</>g(s)) does not depend on 
s if n — > oo (and thus e — > 0): 



Here, A 21 (s) = A 2 (s) - Ai(s) and v £ Q(l/(n + l) 2 ). 

Also by inspection of Eq. (20), A 2 i(s) takes its min- 
imum value at s = 1/2 in the region a < s < 1/2. In 
this case a' = fi' = 2a and the solutions for p at s = 1/2 
satisfy 



lim \<f)q(s)) = 



(26) 



sin((n + 2)p) - 2asin((n + V)p) 
sin((n + l)p) — 2asin(np) 



= ±V2 



(31) 



7 



This is satisfied if p = ir/(n + 1) for a = l/v°j which 
corresponds to the eigenvalue of the second excited state: 



A 2 (l/2) 



1 



(32) 



or, equivalently, x(l/2) = for this solution. The 
eigenvalue Ai(l/2) can be well approximated if we let 
sin((n + 1 + r)p) ps -n(r + x(s))/(n + 1) in Eq. ([3l]) and 
using the minus sign in its rhs. With this approxima- 
tion, the equation becomes linear in x(s) and we obtain 
x(l/2) = -2.82. It means that 



Ai(l/2) 



cos 



1 - 



2.82 



(33) 



with (i e C(l/(n + l) 3 ). 

Therefore, the spectral gap between the second and 
first excited states, for a < s < 1/2, satisfies 



A a i(s) > A a (l/2) - Ai(l/2) 
> c'/n 3 , 



(34) 



for a constant c' « 2.827T. 



